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Abstract 

This talk is based on the Very High Energy Gamma Ray Astronomy 
observations planned to be carried out at Pachmarhi in the central Indian 
state of Madhya Pradesh using the well known atmospheric Cerenkov tech¬ 
nique. The development of a ground based array of 25 Cerenkov telescopes 
is currently underway at Pachmarhi situated at an altitutde of about a kilo¬ 
meter. Using this array it is proposed to sample the Cerenkov light pool 
at various distances from the shower core in order to estimate the lateral 
distribution parameters of the shower. Simulation studies have shown that 
these parametrs would enable one to distinguish gamma ray initiated show¬ 
ers from those by cosmic ray charged particles, thus significantly improving 
the signal to noise ratio. After summarizing the genesis of VHE gamma 
ray astronomy in our institute we will discuss the scientific motivation of 
this concept of enriching the gamma ray signal as compared to the stan¬ 
dard imaging technique. The current status of the detector development 
and the expected results will be presented. 


1 Introduction 

The gamma ray region of the electromagnetic spectrum has been the last to 
be successfully exploited as a channel for astronomical investigation. Although 
it was realized as far back as 1958 (Morrison, 1958), the detection techniques 
were difficult and the fluxes were low. The chief motivation for the study of 
7 -rays is that they are produced in the high energy particle interactions. The 
close association between cosmic rays and y-rays makes y-ray astronomy a special 
branch of cosmic ray astrophysics. A careful study of the intensity spectrum of y- 
rays and their time variability may reveal the physical conditions and acceleration 
mechanisms of cosmic rays with in the source. 

*A talk delivered at the Tata Institute of Fundamental Research during “Perspectives in 
High Energy Astronomy & Astrophysics”, an International Colloquium to commemorate the 
Golden Jubilee year of TIFR (August 12-17, 1996). 
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Before ~1980, “Gamma Ray Astronomy” was interpreted as a branch of Space 
Science and presumed to terminate where satellites ceased to be useful, i.e. ~1 
GeV. Today the situation is quite different. Very high energy gamma ray as¬ 
tronomy, covering the energy range 10^-10^ GeV, has generated lot of interest 
as many new groups all over the world are entering the held. Further the held 
has come of age and one can say with certainty that some VHE 7 -ray sources do 
exist. 


2 Genesis of 7 -Ray Astronomy in TIFR 

At the Tata Institute of Fundamental Research we have a long history in High 
Energy Physics. The interest in air Cerenkov radiation started in the early 60’s 
when we were engaged in the study of characteristics of high energy nuclear 
interactions in the tens of GeV range - using cosmic rays at mountain altitudes 
before the hrst GERN accelerator. The experimental set up at the Gosmic Ray 
Laboratory at Ootacamund in the Southern Indian State of Tamilnadu consisted 
of a triple arrangement comprising of an air Cerenkov counter, a multiplate cloud 
chamber and a total absorption spectrometer one below the other. The hadron 
interaction was recorded in the cloud chamber, its total energy was measured 
by the TASS which was essentially a total ionization calorimeter and the air 
Cerenkov counter above acted as a threshold counter and served the important 
purpose of distinguishing between pions and protons in the energy range 10 - 
40 GeV. The Cerenkov counter was 4 m high and the light was collected by a 
parabolic search light mirror of Im diameter. 

Immediately after the discovery of Pulsars, a series of experiments were de¬ 
signed in 1969 to search for very high energy y-rays from Pulsars using the atmo¬ 
spheric Cerenkov technique (Chatterjee et al, 1970, 71). The first atmospheric 
Cerenkov telescope consisted of two parabolic mirrors and the activity continued 
with a progressively increasing sensitivity. The collection area of mirrors were 
progressively increased to 6.4 m^ (using 10 mirrors of 0.9 m diameter) in 1977 to 
20 m^ in 1979, by adding 8 more mirrors of 1.5 m diameter. The mirrors were 
deployed in the form of a compact array and the signals from neighboring mirrors 
were added linearly to form a bank. A time coincidence between the banks was 
used as the trigger. Later majority logic between different banks was used to 
generate an event trigger. 
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Figure 1: Phasogram of PSR 0531+21 for VHE 7 -rays of energy > 4.5 TeV. Two 
signal peaks separated by the characetristic 0.42 of the pulsar phase could be 
seen. 

3 Past land marks 

3.1 Compact Array: 

3.1.1 Crab Pulsar Observations 

During the 1976-77 observing season the Crab pulsar was observed for 27 hours 
using 10 mirrors of 0.9 m diameter with a mask of 1° FWHM. The event arrival 
times were converted to the arrival times at the solar system barycenter. Using 
the contemporaneous radio ephimeris for Crab, the pulsar phasogram is derived. 
Two peaks of 3.6a and 2.2a signihcance with a separation of 0.42+0.03 as seen 
in the MeV-GeV data were observed. Figure 1 shows the phasogram of Crab 
pulsar in the TeV energy range. The estimated steady flux of VHE 7 -rays of 
energy above 4.5 TeV is estimated to be 2x10“^^ photons cm“^s“^. Ever since 
no observations (before applying any cuts) carried out at Ooty or Pachmarhi 
showed any evidence for pulsed steady emission from this pulsar (Gupta et ai, 
1978). 

3.1.2 Vela Observations 

Vela pulsar (PSR 0833-45) is a gamma ray pulsar which is the brightest 7 -ray 
source in the GeV range. Despite the large angle this source could just be ob¬ 
served from Ooty. This source was observed from Ooty for the hrst time in 
February, 1979. The observations were continued in February-March, 1979. The 
data from 33 hours of observation showed two distinct peaks in the phasogram as 
shown in hgure 2. The integral flux of Gamma rays of energy > 14 TeV from this 
source is estimated to be (5+1) 10“^^ photons cm“^s“F Because of inadequate 
time keeping the absolute pulsar phase could not be established (Bhat et ai, 
1980). 
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Figure 2: Phasogram of about 5000 showers from the direction of Vela Pulsar 
(PSR 0833-45). The data were collected during February-March 1979. 

3.2 Distributed Array: 

The above observations were carried out at the back yard of the Cosmic Ray 
Laboratory where the Cerenkov telescopes were initially deployed in the form of 
a compact array. After a series of tests involving the timing of Cerenkov front 
we showed that it is possible to determine the arrival direction of the shower 
accurately by the well known method of triangulation. This would allow us to 
improve the signal to noise ratio by being able to reject off-axis showers which are 
due to cosmic rays. Hence it was decided to build a distributed array of Cerenkov 
telescopes in order to make use of this idea in 1980. Because of paucity of space 
we had to move and build the array at place about 11 km away. The new array 
consisted of 8 large mirrors at the center of a circle of radius 55 m in the form 
of a compact cluster while 10 small mirrors were situated symmetrically on the 
circumference. While central cluster was used for generating the trigger the outer 
mirrors were used for timing information. The uncertainty in the arrival direction 
estimate is ~ 0.3° if > 6 timing telescopes are triggered in an event. Most of 
events however had only about 3 timing mirrors triggered resulting in the angular 
accuracy of 0.6-0.7° for vertical incidence. One of the important drawbacks of a 
distributed array is that since the mirrors are distributed the energy threshold is 
invariably higher compared to that for a compact array. The remaining 2 small 
mirrors were operated at the old site (Gupta, 1983). 

An added advantage of the distributed array is that the number of triggered 
timing telescopes is an indication of the energy of the primary thus enabling one 
to provide an energy spectrum of gamma ray signal from a celestial source. This 
is a unique capability available for the hrst time in the history of the held of TeV 
7 -ray astronomy. 
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PHASOGRAMS OF VELA 
DATA WITH CUTS 



Figure 3: Phasogram of the 3 datasets for Vela Pulsar after applying cuts to 
reject higher energy showers. 

3.3 Steady signal from Veal Pulsar 

Observations on this source continued almost every year following the initial suc¬ 
cess. The data taken during 1979-80, 1982-83 and 1984-85 with differing exposues 
and energy thresholds (9.3, 4.7 & 5.4 TeV respectively) were independently ana¬ 
lyzed. The phasograms generated using contemporaneous radio ephimeris of this 
pulsar showed a small but persistent peak at the same phase bin over the years 
which increased in strength when a cut is applied to the data to reject 50 % of the 
higher energy events based on the pulse height data (see hgure 3). The combined 
a phasogram showed an impressive 4.0 cr signal whose position coincided with the 
main optical pulse of the pulsar (Bhat et al, 1987). 

Since the different data segments had different energy thresholds one could 
derive an energy spectrum of VHE y-rays for this source for the hrst time. The 
slope of the integral energy spectrum in this energy range is estimated to be 
-(2.5 ±0.3). 

3.4 A GRB from Crab Pulsar? 

On the other hand the Crab pulsar did not show any evidence for steady emission. 
However there was an interesting detection of a transient emission lasting for 
about 15 mins from this source on January 23, 1985 at 1711 hrs UT from this 
source. The phasogram (shown in hgure 4) in the TeV energy range showed a 
strong ( 5.1 a) peak that coincided with the radio main pulse amounting to a 
gamma ray hux of (2.5±0.6)10“^'^ photons cm“^s“^ at E.y > 1.2 TeV (Bhat et al, 
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Figure 4: Phasogram of PSR 0531+21 in TeV 7 -rays during the 15 min 1711-1726 
UT on January 23, 1985. 

1986). This is the second detection of flaring activity of a similar time scale and 
strength in Crab pulsar in the TeV range, the hrst being in 1982 detected by the 
Durham group (Gibson et al, 1982). 

3.5 TeV 7 -rays from Geminga 

Geminga (PSR 0633+23) is the strongest 7 -rays source in the sky in the medium 
energies. A periodicity of 237 ms was seen in the MeV 7 -ray signal from this 
source (Hermsen et al., 1992). This source was observed from Ooty during 1984- 
85 using the distributed array amounting a total of 30 hrs. exposure. When 
searched for a periodic signal in this dataset a modest signal consisting of two 
peaks separated by 0.5 of the pulsar phase was seen, thus establishing the first 
ever detection of TeV 7 -rays from this source. The peak positions coincided with 
those seen at MeV energies. The estimated integral flux of TeV 7 -rays above 0.8 
TeV is (2.1+ 0.8) 10“^^ photons cm“^ s“^(see Vishwanath et al, 1993a). 

4 High Energy Gamma Ray Observatory (HEGRO) 

Since the inception a persistent effort was made to increase the sensitivity of 
the 7 -ray telescopes by (a) improving the hardware deployed in the observations, 
like using better and faster data recording system & increasing the mirror area 
and (b) by developing more techniques to help us reject the background events 
more efficiently like measuring the direction of the primaries as well as measuring 
pulse heights which are proportional to the Cerenkov light produced by them. 

There is yet another parameter which can improve the sensitivity further viz. 
the exposure time. Ooty weather being what it is the total exposure one could 
get on a source was severely limited by poor observing conditions. Hence it was 
soon decided to move the entire array to a new site at Pachmarhi ( longitude: 
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Figure 5: The trigger rate and the chance coincidence rate for every 2 minutes 
for Her X-1. The right hand scale shows the percentage excess over the cosmic 
ray background. 


78° 26' E, latitude: 22° 28' N & altitude: 1075 m) in the central Indian state of 
Madhya Pradesh. After a systematic study of meteorological data of this place 
was chosen which is expected to allow much longer exposures on sources of our 
interest. All the mirrors were deployed in the form of a compact array initially to 
achieve as low an energy threshold as possible for y-rays. The observations were 
resumed within 2 months of shifting. 

4.1 A TeV 7 -ray burst from Her X-1 

The hrst source being observed from this site was the Her X-1 which by then 
has been accepted to be a source of TeV y-rays. During a run on April 1986, a 
large increase in event trigger rate was observed lasting for a period of about 14 
minutes and no such increase was observed in the chance coincidence rate which 
was also being monitored in real time. This is by far the largest busrt (42 a) 
ever seen in this energy range from any source(see hgure 5). During the bursting 
period the time averaged y-ray flux was estimated to be 1.8 xl0“® photons cm“^ 
s“^ for > 0.4 TeV (Vishwanath et ai, 1989). This source has been known 
to be active during certain epochs during its 35 day cycle observed from X-ray 
observations. 
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Figure 6 : 7 -ray phasogram for 428,639 events from the direction of PSR 0355-1-54 
showing a 4.3 a peak at a phase of 0.53, with respect to the radio pulse. 

4.2 Discovery of PSR 0355 + 54 

The motivation to observe this short period (156 ms) radio pulsar was its ex¬ 
tremely low timing noise which might indicate the presence of an outer gap in 
the pulsar magnetoshpere near the light cylinder. The presence of an outer gap 
may mean the existence of energetic particles which in turn could potentially 
produce TeV 7 -rays (Cheng et ah, 1988). In addition, two giant glitches were 
detected from this pulsar in 1986 (Lyne, 1987). From our observations of the 
Vela pulsar which shows frequent glitches, there seems to be an increased 7 -ray 
activity following a glitch. Hence we observed this source using a compact array 
during the end of 1987. During an observation period of ~ 25 hours a convincing 
evidence of a steady emission of TeV 7 -rays was detected from this source for 
the first time. The phasogram, shown in figure 6 , derived using contemporaneous 
pulsar elements showed a 4.3 a peak at a phase of 0.53 with respect to the radio 
pulse amounting to a 7 -ray flux of 7.9 x 10“^^ photons cm“^ s“^ at E.^ > 1.3 TeV 
(Bhat et al, 1990) 

4.3 Distributed Cerenkov Telescopes 

In September 1988, the Pachmarhi array was modihed in order to observe tran¬ 
sient emissions. The available hardware was separated into 4 groups A,B,C & 
D and were deployed at the corners of a rectangle 80 m x 90 m. Each group 
consisted of 2 large (1.5 m diameter) mirrors and 2 smaller (0.9 m diameter) 
parabolic mirrors. Transient emission, if any, could be independently detected by 
these telescopes. Two of them could be monitoring the off-source region during 
each observation. At the center of the rectangle was a fifth array (E) consisting 
of 4 smaller reflectors in order to study the energy threshold dependence of the 
episodal emission. 





CRAB PIII.SAR (2 JAN 19R9) 



Figure 7: Top; Phasogram of PSR 0531+21 during the 5 min burst. Bottom: 
Phasogram for the 1 min interval during the burst. M, I refer to the radio main 
and inter-pulse positions respectively. 

4.3.1 A second Crab burst 

This source was observed continuously during 1988-90 period using this modihed 
array. On January 2, 1989 a 5 minute transient emission was seen from this source 
by the 5 telescopes independently all of which were observing the source. The 
signal strengths were consistent with the different y-ray energy thresholds. This 
transient emission from the Crab pulsar, with a chance probability of 9x10“^, was 
the second observed by our group. The time averaged flux during the 5 minutes 
was 1.4 X 10“^*^ photons cm“^ s“^ for E.^ > 2.8 TeV. The pulsar phasogram 
during this episode showed a strong 6 a peak coinciding with the radio main 
pulse as shown in hgure 7. Some activity in the inter-pulse region also could be 
seen in the phasogram (Acharya et aL, 1992). 

In addition to the TeV y-ray sources mentioned above several other sources 
have also been observed both from Ooty and Pachmarhi which however did not 
show any evidence for the emission of TeV y-rays. They are the nearest pulsar 
PSR 0950+08, the short period X-ray binary 4U0115+63, the millisecond pulsar 
PSR 1957+20 etc. 
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Figure 8 : 


5 Pachmarhi Array of Cerenkov Telescopes 

5.1 Polemics 

In a constant effort to improve the signal to noise ratio further, it was realized 
that one has to reduce the background due to hadronic cosmic ray showers. 
Monte-Carlo simulations have shown that differences between hadron induced 
and 7 -ray initiated cascades could be exploited to reject cosmic ray initiated 
showers (Hillas and Patterson 1987, Rao and Sinha 1988, Knapp and Heck, 1995). 
The simulations have also shown that electromagnetic cascades are flatter in the 
lateral distribution of Cerenkov photons and more compact in the angular size of 
the Cerenkov images than those initiated by cosmic ray primaries. 

Several groups have used some of the above parameters to reject the cosmic ray 
background and enhance the signal (Weekes 1989, Baillon et a/.,1994, Tiimer et 
a/., 1985 and Goret et ah, 1993). The Whipple group has detected steady emission 
of TeV 7 -rays from Crab nebula using Cerenkov imaging technique (Vacanti et 
ah, 1991 and Punch et al.,1992) and have established this source as a ‘standard 
candle’ of TeV 7 —rays. 

However, not enough attention has been given to the lateral distribution as¬ 
pect of the atmospheric Cerenkov radiation to reject cosmic ray initiated back¬ 
ground. As shown in hgure 8 , angular imaging and spatial multiple sampling are 
two complimentary ways to examine the same bundle. However one looks at dif¬ 
ferent types of distinguishing features between 7 -ray and cosmic ray primaries in 
the two techniques. The lateral distribution of Cerenkov photons have a “hump” 
at distances of 120-140 m from the shower core in 7 -ray initiated cascades only 
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(Rao and Sinha 1988). Further, the shower to shower fluctuations in the lat¬ 
eral distribution of Cerenkov photons is much less for a 7 -ray initiated shower 
compared to cosmic ray initiated showers, which are rather ‘bumpy’ due to the 
contribution of Cerenkov light from muons. The average lateral distribution of 
Cerenkov photon densities, as calculated by Hillas and Patterson ( Hillas and 
Patterson 1990), clearly show a flat distribution up to the “hump” region for 
7 -ray initiated showers and a steeper distribution for proton initiated ones. Rao 
and Sinha (Rao and Sinha 1988) have shown that the signal to noise ratio could at 
least be improved twice by exploiting these very differences. Calculations made 
for various array configurations at Pachmarhi showed (Vishwanath et ah,1993) 
that unambiguous identification of 7 -ray initiated showers could be made with 
a large number of detectors sampling the Cerenkov photons. Accordingly, the 
set-up at Pachmarhi is being modified and augmented. 

5.2 The array 

We have decided to pursue the latter technique viz. to measure the lateral dis¬ 
tribution parameters in order to distinguish the photon primaries from the more 
abundant hadronic primaries. The entire available area of 85 m x 100 m at 
HEGRO Pachmarhi has been filled with an array of 25 Cerenkov telescopes. 
Each telescope consists of 7 parabolic (f/d ~ 1) mirrors mounted symmetrically 
making up a total reflector area of 4.4 rn?. These parabolic reflectors are fab¬ 
ricated locally and their optical quality is such that a point source image is < 
1°. The reflectors are equatorially mounted and independently steerable both in 
E-W and N-S directions. The reflector orientation and tracking are controlled by 
a computer automated system (Gothe, et al, 1997) to an accuracy of ±0°.l. A 
fast photomultiplier (EMI 9807B) is mounted at the focal plane of each reflector 
behind a 3° diameter mask. 

The existing large and small mirrors will be deployed in the form of a compact 
array at the center. These will generate an independent trigger and hence will 
work like a compact array within the larger array, having a significantly lower 
energy threshold. From this one will be able to derive an energy spectrum of 
gamma rays from a source 

Event arrival times are derived from a Global Position Satellite (GPS) receiver 
having an absolute time keeping accuracy of ±100 ns. The pulse heights and the 
times of arrival of pulses at each PMT will be recorded using LeGroy ADG and 
TDG modules respectively. A new distributed data acquisition system based on 
several PG 486’s networked together using a LINUX system is currently being 
developed (Bhat, 1996). Event triggers will be generated by each of the 25 
telescopes and tagged for identification. The idea of a distributed data acquisition 
is to minimize the loss of information due to transmission through long coaxial 
cables. Special cables (RG 213) of shorter lengths will be used for analog pulse 
transmission from the phototubes to the nearest signal processing center. Event 
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triggers will be generated from the royal sum of all the analog pulses from all 
the 7 mirrors in a telescope. The measurement of time differences of triggers 
from various telescopes will enable us to measure the arrival direction correct to 
0.2-0.3°. 

Each event data will consist of the event arrival time to an accuracy of ±1 jj,s, 
amplitude and relative time of arrival of pulses at each bank, trigger information 
and other relevant house-keeping informations etc will be recorded at each of 
the 4 held signal processing centers while the relative arrival times between the 
telescopes will be recorded at the central data acquisition system. The data will 
be recorded in the internal hard disks which will be collated off-line. 

A Computer Automated Rate Adjustment and Monitoring System (CARAMS) 
is developed for setting the high voltages on individual phototubes such that a 
preset rate of pulses above the set threshold of 30 mV is generated (Bhat, 1996a). 
This will also monitor the counting rates from all the phototubes through out the 
run. Similarly, the optimum count rate from each photo tube is estimated using 
the count rate variation as a function of the discriminator threshold. A soft¬ 
ware package (Automated Rate Measurement At different Discriminator thresh¬ 
olds, ARMADA) to derive this curve is also developed. This will run prior to 
CARAMS before starting an observation. 

5.3 Photon 8z Energy Thresholds 

If we use a the royal sum of the analog pulses from the 7 mirrors in a telescope 
to generate a trigger, it has been estimated that the trigger rate will be around 
1000/min. The night sky background flux at Pachmarhi was measured to be 
1.2 X 10® photons cm“^ s“^ (Bhat and Mehta, 1996). Using this it is estimated 
that the threshold photon density is ~ 32 7 m~^. Using the estimated photon 
densities through monte carlo simulations at various 7 -ray energies (Ong, 1995), 
this photon density threshold translates to an energy threshold of ~ 350 GeV for 
7 -rays. 

5.4 Technical Progress 

Eight of the 25 telescopes are already in position and the rest will be erected in 
about 3-4 month’s time. 

The new distributed data recording system as well as networking between 
PC’s carrying out signal processing in the held signal processing centers (FSPC) 
are being developed. Each of the FSPC’s will have independent clocks which will 
be synchronized with the GPS through the network. The FSPC’s are constantly 
monitored and controlled from the central station. Also the data recording system 
in the central station communicates with the PC’s which control the high voltages 
on the individual phototubes as well as the PC which controls and monitors the 
movement of the telescopes. 
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We expect to see the first light around October, 1997. 
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